Effect of Deformation on Surface Characteristics of Finite Metallic Crystals* 
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The surface stress and the contact potential differences of elastically deformed faces of Al, Cu, Au, 
Ni, and Ti crystals are calculated within the modified stabilized jellium model using the self - consis- 
tent Kohn-Sham method. The obtained values of the surface stress are in agreement with the results 
of the available first-principal calculations. We find that the work function decreases/increases lin- 
early with elongation/compression of crystals. Our results confirm that the available experimental 
data for the contact potential difference obtained for the deformed surface by the Kelvin method do 
not correspond to the change of the work function but to the change of the surface potential. The 
problem of anisotropy of the work function and ionization potential of finite sample is discussed. 
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I. INTRODUCTION 

The preparation and study of nanometer scale struc- 
tures attracts a considerable current interest for both 
technological and scientific reasons. The free-electron gas 
models are invariably popular tools in the physics of met- 
als P, and low-dimensional structures Q . 

Recently, several authors discussed the definition of the 
surface stress and the chemical potential for the case of 
finite samples 0, H H |L 

In our previous papers [y, Q we proposed a successful 
method to calculate the surface stress and the work func- 
tion of a single elastically deformed finite metal crystal. 
It was applied for simple metal. A fairly tightly bound d 
band that overlaps and hybridizes with a broader nearly 
free - electron sp band characterizes transition metals. 
The cohesive properties of simple, noble and transition 
metals can be calculated from the first principles in the 
context of one-electron picture provided by the density- 
functional theory ll3 • 

The direct measurements, using the Kelvin method 
showed a decreasing/increasing of the contact potential 
difference (CPD) of the elastically tensed/compressed 
metal samples [ill E ES A similar effect on CPD 
was observed at the surface of sample with a nonuniform 
distribution of residual mechanical stress ^5|. The con- 
ventional method of measurement of the work function 
changes versus x-axes strain 0, 0, 0, is based on 
the expression: 



flxx, M Pa 



WKei, eV 



AW, 



K,:l 



W{u xx ) - W(0) = -CPD, 



(1) 



i.e., the work function as if increases for a tensed sample. 
The change AWKei(u xx ) was measured (see FigQJ at the 
side perpendicular y- and z-directions, u xx is the rela- 
tive deformation along x-axes (see FigEJ • These, at first 
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FIG. 1: The strain dependence of external mechanical stress 
Tlxx(u xx ) and work function WKei(uxx), denned by Eq.Q for 
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FIG. 2: Sketch of deformation. 



sight, surprising results mean that the work function in- 
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creases/decreases with uniaxial tension/compression of 
metallic sample. This fact contradicts to another well- 
known observation: the work function of simple metals 
decreases for the transition Al Na => Cs, i.e. with 
decreasing of the electron concentration. 

In the present paper, we report the results of calcula- 
tions of the effect of deformation on the surface energy, 
the work function, and the contact potential difference 
of faces of technologically important metals such as Al, 
Ti, Ni, Cu, and Au using their nominal valence. The 
problem of an accurate definition of the work function is 
discussed. We checked the accuracy of Eq.Q by fully 
self-consistent calculations and showed that it is wrong. 
Eq.Q is incorrect in the diagnostic of strained surface. 



II. NONSTRAINED SURFACE 

Surface energy per unit area and work function are 
the most important characteristics of a metal surface. 
In the framework of density-functional theory the total 
energy of metal is a functional of the nonhomogeneous 
electron concentration n(r), n{r) n in the bulk, h = 
3/47rr 3 , r s o = Z 1//3 ro, r s a and r$ are, respectively, the 
average between electron and between ion distances in 
the nonstrained metal bulk, Z is the metallic valence. 
Using a pseudopotential approach the total energy can 
be written as the sum 

E[n{r)\ =T S + E ex + E cor + E H + E ps + E M , (2) 

where T s is the (non-interacting) electron kinetic energy, 
E ex is the exchange energy, E cor is the correlation en- 
ergy, Eh is the Hartree (electrostatic) energy, E ps is the 
pseudopotential correction, and Em is the Madelung en- 
ergy. The sum of first four terms in Eq. © corresponds to 
the energy of ordinary jellium, Ej. The ordinary jellium 
model provides a useful description for the bulk modulus 
and surface energy of the simple metals only for r S Q 4 
bohr (Na), where bulk jellium is stable. The jellium sur- 
face energy for r S Q < 2 and the bulk modulus for r S Q > 
5 bohr are negative. These deficiencies are removed in 
the stabilized jellium model |l7j . The average energy 
per valence electron in the bulk of stabilized jellium is 
egj = E[n]/N, where TV is a total number of free elec- 
trons, defined by valence and atomic density, 



and the popular expression for the correlation energy [Tsl 



£s.J = £j + w R + £ M , 



(3) 



£ C or[n(r)} 



0.1423 



1 + 0.8293?i(r)- 1 /6 + 0.2068n(r)-V 



I/a- (5) 



Applying the variation principle one can find the Euler 
Lagrange equation for nonhomogeneous electron gas 



(j, = V eff {r) + 



5n(r) ' 



(6) 



where p is the chemical potential of the electrons, and 
effective one - electron potential is 

V eff (r) = p(r) + 5[EeX 5 ^ C ° r] + (SV) face 6(r r>), (7) 

where r' is the radius-vector of the surface, 9(r — r') is the 
step function. The electrostatic potential ip(r) satisfies 
the Poisson equation. 



V 2 ip(r) = -4n[n(r) - p(r)]. 



(8) 



The ionic charge distribution can be modelled by the step 
function p[r) = p0(r — r'), where p = uq/Z. 

The coordinate-independent term (SV)f ace in Eq.© 
represents the difference between the pseudopotential of 
the ion lattice and the electrostatic potential of positive 
background averaged over the Wigner-Seitz cell. This 
term allows one to distinguish different faces of semi- 
infinite samples. The face-dependence of the stabilization 
potential 0,0] reads 



(SV)ws = n Q -—(£ M + w R ) 
cluq 



(9) 



(10) 



where the term 



where do is the spacing between the lattice planes par- 
allel to the surface. From the bulk stability condition 
desj/dho = one can obtained the relation 



(SV)wS — -n -J^(t s + Sex + Ecor), (H) 

ano 

The electronic profile n(r) can be expressed in the terms 
of one-electron wave functions ipi, 



£.J 



t s + £e 



10 



-in 



k F + £ c 



(4) 



consists of the average kinetic, exchange and correlation 
energy per electron, kp — (37r 2 no) 1 / 3 , tq is the radius 
of the Wigner - Seitz cell, and wr = l-Kn^r^ represents 
the average of the repulsive part of the Ashcroft model 
potential, r c is the radius of ionic core, Sm = — 9Z/10ro. 
We employ atomic units throughout (h = m = e = 1) 



N 



(12) 



where the functions ipi satisfy the one-electron wave 
equation 

-\^ 2 Mr) + V e ff ( f ace) (r)ipi{r) = £iipi(r). (13) 
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The set of the Kohn-Sham equations must be solved self- 
consistently. Then the total kinetic energy of electrons 
in EqEJis 



T a [n] 



N 

£< 

i=l 



d 3 rn(r)V ef f(r) 



(14) 



As a rule for semi-infinite metal it is supposed that 
the electronic profile n(r) and effective potential V e //(r) 
vary only in the direction perpendicular to the surface. 
The conventional approach involves introducing periodic 
boundary conditions in the x- and y-directions. Thus, 
only crystal face is in z-direction. In this case the set of 
equations ©, lfT2l and lj*H?|> reduces to 



oo oo 



tp(z) = tp(oo) - 4vr / dz' / dz"[n(z") - p(z")], (15) 



where S is the phase shift of wave function,^ (z) — > 
sinffcz — 5(k)] as z — > — oo, l4//(/ace) is the bulk mag- 
nitude of effective potential, and 



dz[n(z)e xc {n(z)) - ne xc (n)9(-z)], (22) 



-„ = - / dz<p(z)[n(z) - p(n(z))]. (23) 



Putting the electrostatic potential in the vacuum equal 
zero, </?(+oo) =0, one can calculate the work function as 
(e F = fc|/2) 



W z -face = ~M = -Ve/f ~ ^F- 



(24) 



III. DEFORMED SURFACE 



(z) = -j / dk(P F -k 2 ) \Mz)\ 2 



^k{z) = -k 2 ^ k (z). 



(16) 



(17) 



Here the effective potential is taken in a local-density- 
approximation 

V eff [z, n] = 4>(z) + V xc {z) + (6V) face 9{~z), (18) 

where n = n(z). The electrons wave number varies in 
the interval (0, ftp). The solution of the problem can be 
reduced to the iteration procedure by means of the follow 
relation 20] 



,_,!:)= / d z>e~ k ^- z '^{n(z>)-p(z')} + 'fMz), 



where i is the number of iteration. 
The surface energy can be written as 



(19) 



1sj = U + (Vws) J dz[n{z) -n ], (20) 

— oo 

The ordinary jellium components are: 



7s 



1 

2^2 







1 dkk 


J-S(k) 



dzn(z) [V e ff(f ace )(z) - V eff ( face )6(-z)\ , (21) 



The strain dependence of the CPD was measured 
for poly crystalline compressed |ll| and tensed samples 
13, 141. One can assume that a polycrystal is assembled 
from a number of simple crystallites. Thus, qualitatively, 
the problem can be reduced to an analysis of tension or 
compression applied to a single crystal. We consider a 
single crystal in the shape of a parallelepiped with sides 
having equivalent Miller indices. For simplicity, the ma- 
terial of the sample is assumed to have the cubic crystal- 
lographic symmetry. We introduce the coordinate system 
with axes and perpendicular to the sample faces. Their 
areas are equal to and respectively (see FigEJ. We ne- 
glect temperature and dimensional effects that differ from 
the approach used in 0, [2l| . 

Let us first express the average electron density in 
a metal as a function of deformation. For this pur- 
pose, consider an undeformed cubic cell of side length 
on, Og = 47rro/3. In the modified stabilized jellium model 
[a, IS] the metal energy is a function of the electron den- 
sity parameter r su = Z 1 / 3 r 0u , spacing d u between the 
lattice planes perpendicular to z-direction, the Poisson 
coefficient for polycrystal u, the Young modulus Y, and 
deformation u xx . 

For uniaxially deformed cell elongated or compressed 
along the x-axis one can write yj : 



d u =d (l- vu xx ) 



(25) 



where d u is the spacing between the lattice planes per- 
pendicular to the y- or z-directions, do is the interplanar 
spacing in undeformed crystal, 



n [l-(l-2v)u xx ] + 0(u 2 xx ) 



(26) 



is the average electron density in deformed metal bulk, 
v is the Poisson coefficient for polycrystal and the corre- 
sponding density parameter is . 



r s [1 + {l-2v)u xx ] 1/3 . 



(27) 
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The stabilization potential (|llfl must be changed by 



(SV) W s = - 




where 

p = -Qixx + n w + n zz ) = -Yu xx (i - 2v) 

is the external pressure. 

An applied force stimulates the change of a volume 
and reticular electron density at particular faces of neu- 
tral metallic crystal and results in a difference between 
electrostatic potentials of these faces. 

The corresponding anisotropic three-dimensional elec- 
tric field arises due to a transfer of electrons between 
faces. Thus, finite sample sizes lead, in the first place, 
to the fact that the surface tends to equipotentiality. In 
general case, this transfer of electrons does not give any 
possibility to calculate the work function using conven- 
tional approach for semi-infinite sample. However, in the 
special case of the largest face of the metallic crystal, 

A z > A x ,A y , (29) 

the use of the conventional approach (see Eq. I|24|0 in 
fact is adequate Q , A z is the area of upper side at Fig. 
2. As a result, the calculated work function exclusively 
for z-direction corresponds to that for the whole crys- 
tal. The face perpendicular to this direction is not ap- 
preciably perturbed by the transferred electronic charge 
independently from the crystallographic orientation [8j- 
The conventional method involves introducing periodic 
boundary conditions in the x- and y-directions. In this 
case the surface potential depends strongly from the 
atomic packing density. Thus, the only crystal face is 
in -direction. Putting the three-dimensional electrostatic 
potential in the vacuum equal zero, under condition l|29(l 
one can use definition 1)24(1 . where the effective potential 
in the bulk of semi-infinite metal yields the total sample- 
vacuum barrier. Solving the set of Eqs. (|I5|I - (|I8|I and 
using relations (|2>J) ~ ((2EJ>, one can calculate the strain 
dependencies of the surface energy and work function. 

We calculate the diagonal component of surface stress 
forgiven largest face using the following expression 0, 0, 

aif 

t xx = 7 + — , (30) 

where 7 is the surface energy per unit area of the largest 
face of sample. 



IV. RESULTS AND DISCUSSIONS 

We performed our calculations for the work function 
and the surface energy for the case of absence of the 
strain state of metals and for the strain dependences 



TABLE I: Calculated bulk modulus B, surface energy "/face, 
and the work function Wface- For univalent Au (Z = l)the 
results are placed in brackets. The experimental values of B 
(are given), a, Wface for polycrystaline metal, and the Poisson 
ratio v, the Young modulus Y necessary for calculation the 
strain dependences are taken from books |23l I24I |25|. 





B 


Face 


& f ace 


Wface 


Y 


V 




[Mbar] 




[erg/cm 2 ] 


[eV] 


[GPa] 








(100) 


1087 


3.806 


62.5 




Al 


f.565/0.722 


(no) 


1683 


3.643 


71.4 


0.34 






(111) 


939 


4.119 


75.1 








(100) 


1069(395) 


3.792(3.318) 


43.5 




Au 


L5f9(0.202)/0.722 


(110) 


1652(440) 


3.630(3.148) 


81.3 


0.42 






(111) 


924(383) 


4.105(3.478) 


115.0 








(100) 


1376 


4.010 


138.0 




Ni 


2.567/L860 


(110) 


2224 


3.858 


215.0 


0.32 






(111) 


1162 


4.325 


262.0 








(100) 


979 


3.855 


65.8 




Cu 


I.lf3/L370 


(110) 


1295 


3.647 


131.0 


0.35 






(111) 


899 


4.123 


194.0 








(0001) 


1081 


4.205 


145.0 








(100) 


1355 


3.865 


96.1 




Ti 


f.565/0.722 


(110) 


2456 


3.774 


96.1 


0.30 






(111) 


1081 


4.205 


27.8 





Wface(u X x) and jface(u X x) within the range of defor- 
mation -0.01 < u xx < +0.01 for Ni and -0.03 < 
u xx < +0.03 for Al, Au, Cu, and Ti, respectively. The 
positive/negative deformation is equivalent to the ten- 
sion/compression of the largest side of the sample, i.e., 
the decrease/increase of the atomic packing density at 
this side, and the decrease/increase of the mean electron 
concentration and the interplanar spacing in direction. 
Upper side of sample in Figgis suggested as having in- 
dexes (100), (110), (111) or (0001). 

The tends in bulk and surface characteristics of Al, Au, 
Cu, Ni, and Ti have been reproduced from the theory of 
uniform electron gas with the average electron density 
corresponding to the nominal valence by the volume per 
atom. The use of the integer and fractional magnitude of 
a valency leads to a successful calculation of bulk mod- 
ulus B • The concept of the metallic valency can be 
defined to treat simple and transition metals in terms of 
the uniform electron gas. 

The results of our calculations for free and deformed 
faces are summarized in TableQJand Figs. 2 and 3, respec- 
tively. For comparison we perform calculations of bulk 
modulus B = n^d 2 {nQ£sj)/dfiQ, 7 and W for Au (Z = 
1). The data of Table [I] demonstrate intricate picture 
for common description of cohesive and surface metal- 
lic properties in the simple approach. In fact the use of 
nominal valence gives satisfactory results for the surface 
characteristics. We found that the deformation changes 
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dg/Uxx, erg/cm 



Uxx<0 



Uxx>0 



Uxx<0 



Uxx>0 




Al Ni Cu Au Ti 



Ti Au Cu Ni Al 



FIG. 3: Calculated derivative dy/du^x for estimation of the 
surface stress, Eg. 12911 . The left and right parts of the figure 
correspond to the compression (u xx < 0) and the tension 
(u xx > 0) of the sample, respectively: □ - fee (100), A - fee 
(110), • - fee (111), - hep (0001). 



of the work function and the surface energy remain linear 
with respect to the deformation. The strain derivative 
of the work function and the surface energy is positive. 
The values of surface stress component vary appreciably 
within the interval (1.15, 1.75)7/ Qce . Similar results for 
the surface stress (1250 and 1440 erg/cm 2 for Al (111)) 
are yielded by the ab initio [2(| and atomistic [23 cal- 
culations. In contrast to the present work (FigEjJ, the 
strain derivative obtained for Au in J28| was larger than 
the surface energy. 

The relative change of the work function equals ap- 
proximately 1% for the maximal strains (compare Table 
13 and FigBJ. For the case of the compression (u xx < 0), 
the tail of electron profile and, therefore, of the effective 
potential grows steeper in vacuum. For the case of the 
tensile strain (u xx > 0), these coordinate dependences 
have the opposite tendency. A total decrease/increase of 
the work function is determined by a positive/negative 
shift of the effective potential versus strain in the bulk of 
metal (neglecting the deformation dependence Sf(u xx ), 
one can certainly put the shift AW ~ — AV e ff). Our 
calculation mimics the usual work function dependence 
from electronic concentration (for transition Al ==> Na 
=> Cs). 

However, relation Q gives incorrect dependence 
W(u xx ). These results, at first sight, contradict to ex- 
periments m m m 03 . 

We suggest that, qualitatively, the problem is reduced 
to the effect of strain on a single crystal. The experimen- 
tal observations can be explained based on the change of 
the effective potential at the position of the image plane 
z = zq |8(. For simplicity we take z = 1 bohr for faces 
and use 



CPD = AVeff (z ,u xx ). 

At present, we calculated AW and CPD without use of 
Eq.Q. Our Fig0]demonstrates, on the one hand, a good 
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FIG. 4: Calculated change of the work function and the 
contact potential difference for maximal strain of elastically 
deformed largest faces (Fig|5J. The experimental values of 
—CPD taken from fTHl for compressed (u xx = —0.03) poly- 
crystalline Al, Cu, and Au samples, and from |l3L Il4| for 
tensed Al (u xx = +0.03) and Ni (u xx = +0.01) samples (see 
also Table 1 in 29]). The values of for negative deforma- 
tion are extracted from the lowest experimental shifts in con- 
tact potential —CPD (measured in units of /iVxcm 2 /kg) of 
H3 multiplied by Y\u xx \, where Y = §(Yboi + Yuo + Yin). 
Note that the values of and were erroneously equated in 
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qualitative agreement of the calculated values of with the 
experimental data, on the other hand, the inverse depen- 
dence than it follows from Eq. 

Let us introduce £ the ratio of the effective potential 
differences between strained and strain-free samples at 
the surface and in the bulk, 



AVeff {z ,u xx ) CPD 



AVeff (U M 



-AW 



where V e ff{u xx ) = V e ff(z = —oo,u xx ). The calculations 
show that the value £ varies within the interval (-3, -1) 
for the metal faces under study. For Al £ ~ —2.8!!! 

The analysis of experimental data provided the evi- 
dence that it is not adequate to use the Kelvin method 
for the measurement, e.g., of temperature dependence of 
the work function (see Fig. 6 of |301|). 

Our results show: 

(i) the strain changes of the effective potential in the 
bulk and at the surface have the opposite signs, 

(ii) the sign of the deformation effect is independent on 
the material and of the crystallographic orientation. 

In the conclusion let us remember that the density- 
functional calculations of this quantity deal with the 
nonrealistic semi-infinite systems. However, in all ex- 
periments finite samples are used, in general, of arbi- 
trary sizes and shapes. In this context, an important 
remark appears about the conventional definition of face- 
dependent work function. The ionization potential (IP) 
of a finite sample (as a giant molecule) is given by 



IP = E 



N-l 



E 



N 



w 



e 

2C 7 ' 



(31) 
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where En is the total energy of neutral metallic sample 
containing TV electrons. Here W is introduced as a work 
function of the whole sample and C is a capacitance of 
the sample, which is a scalar quantity. It is clearly seen 
from Eq. 131(1 that the work function is not an anisotropic 
characteristic but a scalar quantity, and IP — ► W at 
C — > co. It may be demonstrated by formation of finite 
sample adding consistently atoms each to other. The 
ionization potential of single atom, dimer, etc. (up to 
work function) corresponds to the ionization process of 
each stage of nucleation (up to the solid) [3l| . Thus, the 
conventional definition (|31[) leads to the conclusion that 
the concept of an anisotropy of the work function of metal 
is spurious in principle. Hence, the realistic geometry 
considerations presented in the flri Is^ ] are adequate in 



the case when all sample faces is posses the same atomic 
packing density. 

Our theory can be applied for the explanation of new 
pressure effects in a single-electron transistor 33] . 
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